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, ( Abstract 

We investigate mass transfer during absorption of atmospheric trace soluble gases by a single droplet whose size is 
comparable to the molecular mean free path in air at normal conditions. It is assumed that the trace reactant diffuses 
£^\to the droplet surface and then reacts with the substances inside the droplet according to the first order rate law. Our 
fc^H analysis applies a flux-matching theory of transport processes in gases and assumes constant thermophysical properties 
of the gases and liquids. We derive an integral equation of Volterra type for the transient molecular flux density to a 

00 liquid droplet and solve it numerically. Numerical calculations are performed for absorption of sulfur dioxide (SO2), 
dinitrogen trioxide (N2O3) and chlorine (CI2) by liquid nanoaerosols accompanied by chemical dissociation reaction. It 

' (-H ' is shown that during gas absorption by nanoaerosols the kinetic effects play significant role, and neglecting kinetic effects 
leads to significant overestimation of the soluble gas flux into a droplet during all the period of gas absorption. 
1 < 
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y 1. Introduction 

£>->, Atmospheric aerosols are directly emitted into the at- 
mosphere from natural or anthropogenic sources or can 

1 ^^ Ibe formed in the atmosphere through nucleation of gas- 

phase species. Aerosol nucleation events produce a large 
fraction of atmospheric aerosols. New particle formation 
occurs in two distinct stages, i.e., nucleation to form a crit- 
ical nucleus and subsequent growth of the critical nucleus 
to a larger size (>2 - 3 nm) that competes with capture 
H and removal of the freshly nucleated nanoparticles by co - 



The sources of these gases in the atm osph ere are briefly 
reviewed in ISeinfeld fc Pandid (|2006l ) and lHavden et al 



> 



(2008|. Soluble gas absorption by no ncirculating droplets 
was i nvestigated experimentally by iTaniguchi fc Asano 
( 19921 ) where conditions of noncirculation for falling liq- 
uid droplets were determined by employing water droplets 
with Sauter mean diameter equal to 0.185, 0.148 and 0.137 



agulation with pre-existing aerosols ([Zhang et all 120121 ). 



in . 

f*^ In the continental boundary layer, there are frequent ob- 
Csl servations of the formation of ultrafine aerosol particles 



Gas absorption by stagnant liquid droplets in the pres- 
ence of inert admixtures when both phase s affect mass 
transfer was analyzed by IClift et all (|1978l ). pp. 54-55. 
Scavenging of s oluble gases by single evap orating droplets 
was studied bv lElperin et all (|2007t l2008f ). 



accom panied by the subsequent growth (jKulmala et al 



> l2004h . Gas absorption of soluble trace atmospheric gases 
[by the liquid atmospheric aerosol particles including ul- 
rN trafine particles plays an important role in climate and 
atmospheric chemistry. 

The consequence for the aerosol climate forcing is 
that the cooling can be intensified with increasing at- 
mospheric amount of water-soluble trace gases such as 
HNO3, count eracting the warmin g effect of the green- 



house gases ( Kramer et al. . [2000). Scavenging of 



at- 



mospheric gaseous pollutants by cl oud droplets is a re- 
sult of gas absorption mechanism ([Pruppacher fc Klett 
1997 : Flossmannl 1990h . Gas scavenging of highly solu- 



All aforementioned works considered the trace gas ab- 
sorption in the continuous limit, where Fick's law re- 
lates the flux and the concentration gradient of reac- 
tant. Recently some studies attempt to describe the 
reactant transport in the gaseous phase in the free- 
molecular and the transition regimes, where the droplet 
size is less than or comparable to the mean free path 
of the molecules in the gaseous phase. Comprehensive 
reviews of the results a nd approaches ca n be found in 
Seinfeld fc Pandisl (120061). IClementl (|2007t >. IPoschl et al' 
(|2007l) and iLushnikovl (120121) . The existing attempts 



ble gases by atmospheric water droplets includes absorp- 
tion of HNO3, H2O2, H2SO4, HCL and some other gases. 
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to describe the transitional and free-molecule regimes 
encounter some difficulties in formulating the boundary 
conditions to the transport equations. For example, 
Vesala et al.l ( 200lh used the diffusion equations and a mi- 
croscopic boundary condition for describing the reactant 
transport toward the particle surface. The rigorous ap- 
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proach requires solution of the full transport problem in- 
cluding the reactant transport in the gaseous phase and 
gas-liquid interface by solving the Boltzmann kinetic equa- 
tion and the diffusion-reaction equation inside the droplet. 
To the best of our knowledge this approach was not applied 
as yet. 

In this study we employ a different approach for de- 
scribing gas uptake b y nanodroplets. We modify th e flux- 
matching approach of iLushnikov fc Kulmala (2004) by in- 
cluding the in-particle chemical transformations of the re- 
actant molecules. The results are applied for considering 
the trace gases scavenging by atmospheric aerosols. 

Let us assume that the reactant molecules (A- 
molecules) move toward the droplet which captures them. 
The further fate of reactant molecules depends on the re- 
sults of chemical reactions inside the particle. Let us de- 
note by n± the concentrations of ^4-molecules right out- 
side (n_|_) or right inside the particle surface. Clearly 
these concentrations n± depend on the nature of physico- 
chemical processes at the surface and inside the particle. 
Let be the number density of A-molecules far from the 
particle. It is commonly accepted that the concentration 
difference — n + drives a flux of ^4-molecules towards 
the particle surface. The mass of the particle begins to 
increase and its chemical composition changes. The rate 
of change of the number of A-molcculcs inside the parti- 
cle is equal to the total flux J, i.e. the total number of 
molecules deposited per unit time at the particle surface 
minus the rate of consumption of molecules A by chemical 
reaction inside the particle. Some fraction of A-molecules 
is assumed to escape from the particles. In steady-state 
conditions the flux J can be written as 



J = a(a)(noo - n + ) 



(1) 



where a(a) is the capture efficiency and a is the particle 
radius. Clearly, capture efficiency a depends on the mass 
accommodation coefficient S p . The latter is defined as the 
probability for an ^4-molecule to stick to the particle after 
a single collision. Since the interface and in-particle pro- 
cesses dete rmine the value o f n+ Eq. (TTJ can be rewritten 
as follows ([Lushnikovl . 120121 ): 



,/ = 



a(a)n c 



1 + a(a)ip(a, J) 



(2) 



where ip{a, J) is a function depending on the nature of the 
chemical reaction. In the case of the first-order chemical 
reaction the function tp is independent of J. An example 
of such function is given below. If, however, the chemical 
reaction inside the particle is of the higher order, then J 
is a solution of the transcendent equation ^ . Note that 
we did not yet specify the functions a(a), ip{a, J), and Eq. 
@ contains only rioo. The information on the processes 
at the surface and inside the particle is carried by the 
function ip(a,J). Hence, Eq. ([2]) is quite general. All 
further approximations concern the values of the uptake 
efficiency a(a) and the reaction function ip(a). 



2. Description of the model 

2.1. Preliminary remarks 

The characteristic times of changes of the number den- 
sity of reactant A in gaseous and liquid phases differ by 
several orders of magnitude. In particular, the relaxation 
time inside the micrometer-sized droplet is t l oc a 2 /D L ~ 
10 -3 s, where D L ~ 10 -5 cm 2 /s denotes the diffusivity 
of the reactant molecules in the liquid phase. The re- 
laxation time in the gaseous phase can be estimated as 
t g (x a/v T ~ 10 -8 s, where v T is the mean thermal veloc- 
ity of the reactant molecules, v T ~ I0 2 m/s. Here one can 
use the free molecular estimate because the droplet size is 
of the ord er of the mean free path o f the reactant molecule. 
Following [Seinfeld fc Pandis (2006) the characteristic time 
of diffusion in a gaseous phase, corresponding to the time 
required by gas-phase diffusion to establish a steady-state 
profile around a particle, can be alternatively estimated 
as t g cx a 2 /D G ~ 10~ 8 s. As can be seen from these es- 
timates, the characteristic time of diffusion in a gaseous 
phase, t g , is much smaller than the characteristic time 
of diffusion in the liquid phase, t l , which is required for 
a saturation of the droplet by soluble gas (i.e., r G <C t l ). 
Therefore for the large values of t (t ^> r G ) it is reasonable 
to assume that concentration profile in a gaseous phase 
in the transitional r egime and the flux attain their quasi- 
steady state values (jWagnerl 1198211 and are determined by 
Eq. <jTJ) (or in more general form by Eq. @). 

The molecular mean free path in air at normal condi- 
tions is £ ~ 65 nm , i.e., it is comparable to the sizes of the 
nanometer droplets. This implies that the motion of the 
reactant molecules cannot be described as the Fickian dif- 
fusion, and one must apply the Boltzmann kinetic theory. 
However, solving the Boltzmann equation analytically or 
numerically is a formidable task. 

The idea to apply the flux-ma tching appro ach in the 
aerosol kinetics was pioneered by iFuchsl (1964). His rea- 
soning was quite simple. At large distance from the droplet 
the reactant transport can be described by the diffusion 
equation. In the vicinity of the droplet at the distances 
of the order of I or less the collisions with the carrier gas 
do not hinder the reactant transport. Consequently, inside 
the region a<r<Roc£(R is referred to as the radius 
of the limiting sphere) the reactant molecules move in the 
free molecule regime. T he value of R mu s t be f ound from 
different consideration. iFuchs fc Sutugm] fll97lh proposed 



to determine this value f rom the num erical solution of the 
BGK equation (see e.g., Sahni, 19661) . An improved ver- 
sion of the Fuchs interpolat ion formula (jFuch s. 1964) was 
obtained by lLovalkal (Il982h in near continuum regime by 
solving Boltzmann equation using momentum method. 



2.2. Trapping el 
The latest 



jciency 
modification 



of the 



(jLushnikov fc Kulmalal . l2004t ILushn ikov 



Fuchs theory 

20121 ) in- 



cludes the solution of the diffusion equation with a fixed 
flux J in the diffusion zone r > R, the solution of the 



2 



collisionless Boltzmann equation in the free molecular 
zone r < R, and determining the radius of the limiting 
sphere from the condition of equality of the fluxes in 
both zones. The expressi on for a(a) was found by 



Lushnikov fc Kulmalal ([2004) for n + = and S p = 1 



a(a) 



2tTQ Vn 




(3) 



where D G is the reactant diffusivity in the carrier gas. The 
extension of this for mula to the case n+ > and S p < 1 
reads (for details see Lushnikovl 120 12h : 



a(a) — 



•} p na vt 



1 + ^ 




(4) 



The radius R of the limiting sphere is found from the 
condition of the equality of flux in the diffusion region 
to the flux from in free molecular regio n. This c ondition 



yield the radius of the limit ing sphere (jLushnikovl [2012 
Lushnikov fc Kulmalal . 120041 ): 



R 




(5) 



It must be noted that R is independent of S p and n + . 
The spherical surface with the radius R separates between 
the zones of the free-molecular and the continuous flow 
regimes. The value of 2D a /v T is of the order of I, the 
reactant molecule mean free path. Hence, if a ~ t or less, 
then the radius R is of the order of R ~ I. 

The concentration profile of the reactant n(r) in the 
gaseous phase inside the limiting sphere a < r < R is 
continuous at r = R together with its first derivative and 
is giv en by the following formula (jLushnikov fc Kulmalal . 
2004h : 

a(a) \ b(r) 



n{r) 



n x - n + 



1 



4ttD g R J b(R) 



where 



b(r) 



1-Ml 



-5 



Outside the limiting sphere at r > R 

n(r) — n + a{a) 



1 



AttDt 



(0) 



(7) 



(8) 



Note that the number density n(a) is always larger than 
n + . The form ula for the concent ration jump at the particle 
surface reads (|Lushnikovl 120121 ) : 



A Q = n(a) - n + = (n^ - n + ) 1 - 



a(a) \ b(a) 



AtiD g R J b(R) 



(9) 



Inspection of Eq. (|4j) shows that when av T /2D G ^> 1, 
a (a) = 47raZ? G and Eq. ^ recovers the Maxwell's equa- 
tion for the molecular flux in the continuum regime J c : 



J c = 4:iraD G (n c 



(10) 



In Fig. [T] we showed the dependence of the ratio J/ J c 
vs. Knudsen number Kn for differen values of the ac- 
commodation coefficient S p (S p was assumed to be 0.1, 
0.2, 0.5 and 1.0 and D G ~ 10 -5 m 2 /s). As can be seen 
from this plot the role of the kinetic effects can be sig- 
nificant for the Kn > 0.1. Comparison of mass trans- 
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Figure 1: Ratio of the molecular flux to the molecular flux in a 
continuum regime, J/ J c , as a function of Knudsen number Kn. 



fer r ates as funct i on of Kn predic t ed by different theo- 
ries dFuchsl. Il964t iFuchs fc Sutueinl . Il97ll: lLovalkal Il982 ; 
Lushnikov fc Kulmalal . I2004T) are shown in Fig. [2l As can 
be seen from Fig. [5] all approaches yield approximately 
the same results for small Kn numbers, Kn ^0.1, and for 
Kn > 10. 



2.3. Inner diffusion-reaction equation 

Consider now the effect of the first order chemical reac- 
tion, e.g., chemical reaction dissociation, inside the droplet 
on the reactant flux towards the droplet. Neglecting re- 
combination the number density of reactant molecules in- 
side the particle ni(r, f) is governed by the linear diffusion- 
reaction equation: 



dn L 
dt 



D L An L - \n L , 



(11) 



where til = til (r, t) is the number density of reactant 
molecules inside the particle, Dl is the reactant diffusivity 
in the liquid phase, and A is the dissociation rate. Equation 
(jlip must be supplemented with the initial condition, 



n L (r,0)=0 



(12) 



0.8 

-3° 0.6 
0.4 
0.2 



0.001 0.01 0.1 1 

Kn 




Loyalka(1983) 
■ Lushnikov&Kulmala (2004) 
- Fuchs (1 964) 

Fuchs&Sutugin (1971) 



100 



Figure 2: Ratio of the molecular flux to the molecular flux in a con- 
t inuum r egime , J/ J c , as a function of Knudsen n umber K n: 
iFuchj dl964T ); I Fuchs fe Sutuginl 1)197111 ; - ■ - - lLovalkal 



(I1982T) ; - •• - - iLushnikov fc Kulmalal ^OmT Accomodation coeffi- 



(no reactant inside the droplet at t = 0)and the boundary 
conditions: 

dn L (r,t) 



and 



Or 



D, 







dn L (r, t) 



Or 



(13) 



(14) 



where j = j(t) is the flux density j — J/(4ira 2 ), J is the 
total flux. 

Equation (jlip with the initial and boundary conditions 
(|12p - (H3J) can be solved by the method of separation of 
variables (see Appendix A) . The result is the concentration 
profile riL(r,t) as a linear functional of j(t). Substituting 
the determined concentration distribution nj, (r, t) into the 
boundary conditi on. Eq. (1141) . yields the integral equation 
of Volterra type ilAuelblatl . bilO-s; for j(t): 



a(a) 
Ana 2 



H J S(t-t')j(t')dt' 
o 



where 



(15) 



(16) 



n> 



and [i n is the infinite set of the roots of the following tran- 
scendental equation: 



[i — tan(/i) . 



(17) 



Equation ()15|) can be rewritten in the dimensionless form 
for the dimensionless flux j*(r) ~ j(t)47ra 2 /a(a)n 00 , r = 
D L t/a 2 : 

T 

fir) = 1 - g(a) J S*(r - r')f(r')dr', (18) 
o 



3a{a)H 

where g(a) = and 

A-kciDl 



f \ — „-Di>(t-t') 



S*(t — t ) = e 



71—1 



(19) 



In Eqs. (TT8"]) and (fl^)) Da, = Xa 2 /Dl is Damkohler number, 
T~L = (Ha'TIT)~ 1 and Ha is Henry's law constant, 1Z is 
universal gas constant and T is the temperature in the 
gaseous phase. 



3. Numerical method 

For the solution of the integral equation Eq. (jT5J) we use 
the method based on the approximation of the integral in 
Eq. (|18l) using some quadrature formula: 

/ F(x)dx = ^2A j F(x j )dx + R n (F), (20) 
j=i 

where € [a, b], j = l,2,...,n, Aj are the coefficients 
associated with a family of quadrature rules and R n (F) is 
a corresponding residuum. Taking successively x = Xi (i = 
1, n) and using the quadrature formula after discarding 
the terms R n (Fi) (i — l,...,n) we obtain the following 
system of linear algebraic equations: 



3i ~ 3(a) E A i S *ijfi = 1 (< = !» -> «)■ ( 21 ) 

The solution of Eqs. (|2ip yields the approximative value 
of the unknown function j$ at the mesh point Tj. The 
system of Eqs. (|2T|) can be written in the following form: 



- £ .l ; A,,./; + (1 - i* = 1, (22) 

i=i 

where = g(a)S*j. In a matrix form the system of Eqs. 
(|2"l"T) can be written as follows: 



/ 1 - AiJTn 

Alif 2 i 1 - ,42^22 

/ 1 \ 

1 

V i J 



( x \ 

3*2 



1 - A n K nn j \ j* j 



(23) 



Using the unequally spaced mesh with an increment 
hi = Ti — 7"j_i, i = 2, and applying the trapezoidal 



integration rule Eqs. (|21[) yield the following recurrence 
equations: 



31 



32 = 



1 - ^K 22 
2 



3=2 



Kijj 



(< 



1 



(24) 



Equations (|24p are valid in the case when hi ^ — . 

In the numerical calculations we spaced the mesh points 
adaptively using the following formula: 



n = ti+(t n - t\) 



1 — cos 



7T i — 1 

2 iV- 1 



( ? ; = i,2,...,iv). 

(25) 



In Eq. (|2"5"j) A" is the chosen number of mesh points, t\ and 
r^v are the locations of left and right boundaries of time 
interval, respectively. 

4. Results and discussion 

Using the suggested model the calculations were per- 
formed for sulfur dioxide (SO2), dinitrogen trioxide 
(N2O3) and chlorine (CI2) absorption by water aerosol 
particles. In order to validate our model we compared 
the results obtained using the suggested model with the 
results obtained in o ur previous study fo r large droplets 



(Kn < l)(see e.g., lElperin et all [2008). The calcula- 



tions were performed for the SO2 absorption by a non- 
evaporating water droplet of 10 fim in radius. The con- 
centration of sulfur dioxide in ambient air was assumed to 
be 0.01 ppm. The calculations showed that the time of the 
complete saturation of droplet by sulfur dioxide estimated 
using the suggested model is 0.08 s, while the time of 
complete saturation of droplet by sulfur dioxide estimated 
using our previous model is w 0.1 s. These calculations 
demonstrate that the results obtained by both models are 
in fairly good agreement. 

The results of calculation of the total mass flux of sulfur 
dioxide as a function of time are shown in Figs. [3] and HI 
The calculations were performed for various radii of water 
aerosol particle (from 0.5 /im to 1.0 fxm, 0.07 < Kn < 0.14, 
see Fig. Hand from 50 nm to 100 nm, 0.7 < Kn < 1.42, 
see Fig. gj). 

As can be seen from these plots for small and moderate 
size droplets the flux of absorbate decreases rapidly at the 
initial stage of gas absorption and approaches asymptoti- 
cally zero at the final stage of the process. The vanishing 
flux of the absorbate implies the stage of saturation of 
droplet by gas. 
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Figure 3: Total molecular flux of sulfur dioxide as a function of time 
(droplet radii a = 1.0, 0.8 and 0.5 fim) 
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Figure 4: Total molecular flux of sulfur dioxide as a function of time 
(droplet radii a = 100.0, 80.0 and 50.0 nm) 



The results of calculation of the dimensionless flux 
J*( T ) — j (t)4:TT a 2 /a(a)n oo as a function of dimensionless 
time r = D L t/a 2 for different gases such as sulfur diox- 
ide, chlorine and dinitrogen trioxide are shown in Fig. [SJ 
Larger values of mass flux at the later stages of gas ab- 
sorption for N2O3 in comparison with SO2 and CI2 ab- 
sorption can be explained by large values of the constant 
of chemical reaction for N2O3 gas absorption in water 
(An 2 o 3 = 1.2 • 10 4 s"\ A so , = 10~ 3 s-\ and A C i 2 = 13.3 

Intensive depletion of the dissolved gaseous N2O3 in a 
water droplet due to chemical reaction leads to the de- 
crease of N2O3 concentration in the bulk of a water droplet 
and to increase of the concentration gradient at the inter- 
face in a liquid phase. These both factors increase mass 
transfer coefficient in a droplet and increase the driving 



5 
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Figure 5: Dimensionless molecular flux density j*(r) as a function 
of dimensionless time r (radius of the droplet R = 100 nm). 



force of mass transfer in liquid. 

As it was mentioned above in the case of large droplets 
(Kn <C 1) the capture efficiency (see Eqs §5§ and can 
be expressed by a(a) — AitaDc- Consequently the flux 
of soluble gas in a gaseous phase is expressed by equation 
similar to Maxwell's equation (see Eq. (ITU1) '). In Fig. [5] we 
showed the results of calculation of the dimensionless flux 
j*(r) of sulfur dioxide as a function of the dimensionless 
time r. It was assumed that concentration of sulfur diox- 
ide in a gaseous phase is equal to 1 ppb, temperature in 
a gaseous phase 298 K and radius of water nanoparticle 
is equal to 10 nm. The dashed line presents the results 
of calculation when the capture efficiency was calculated 
using Eq. Qj. In our calculations we assumed that the 
accommodation coefficient S„ = 1. 




X 



Figure 6: Dimensionless molecular flux density j*(r) as a function 
of dimensionless time r (plane line - a(a) = iiraDc, dashed line — 
a(a) calculated using Eq. ((3j ) . 



The solid line present the results of calculation without 
using the kinetic approach, and thereby the capture effi- 
ciency was assumed to be equal a(a) — AttciDq. As can be 
seen from these plots neglecting kinetic effects in the case 
of gas absorption by nanoaerosols can lead to the essential 
overestimation of mass flux. 

Dependence of the average concentration of soluble sul- 
fur dioxide in a droplet vs. time is shown in Figs. [7] and 

m 

1.4 lO" 6 | 1 
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Figure 7: Concentration of the dissolved SO2 in the bulk of a droplet 
as a function of time (radii of a droplet 0.5, 0.8 and 1.0 fim). 
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Figure 8: Concentration of the dissolved SO2 in the bulk of a droplet 
as a function of time (radii of a droplet 5.0, 10.0 and 20.0 nm). 

Calculations were performed for water droplets with the 
radii 0.5, 0.8 and 1 /jm (Fig. [7J and for water droplets of 
radii 5, 10 and 20 nm (Fig. [5]). In these calculations we 
employed the kinetic approach by using Eq. ^ for the 
capture efficiency a{a) with S p = 1. In Fig. [9] the depen- 
dence of average concentration of SO2 in a droplet vs. time 
was calculated using the kinetic approach (solid lines) and 



G 



neglecting kinetic effects (dashed lines). Calculations were 




' ~~ a(a) - 4 71 a D G 
— a(a)-Eq. (3) 
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Figure 9: Effect of Knudsen layer on temporal evolution of concen- 
tration of the dissolved SO2 for droplets with the radii 5.0 and 10.0 
nm. 



performed for the droplets with the radii 10 and 5 nm. 
As can be seen from these plots neglecting kinetic effects 
leads to the significant overestimation of concentration of 
the dissolved gas in a droplet during all the period of gas 
absorption. Clearly, when the duration of gas absorption 
t — > 00, both approaches yield the same result for the 
magnitude of the dissolved gas concentration. 



5. Conclusions 

In this study we developed a model for absorption of 
soluble trace gases by nanoaerosols taking into account 
dissociation reaction of the first order in a liquid phase. 
In the case when radius of the particle is comparable with 
the mean free path transport of reactant molecules cannot 
be described by Fickian diffusion. However, application 
of the flux-matching theory allowed using transient diffu- 
sion equation with the kinetic boundary conditions for the 
description of gas absorption by nanoaerosols. Transient 
diffusion equation was solved analytically by the method 
of separation of variables. We derived linear integral equa- 
tion of Volterra type for the transient mass flux to a liquid 
droplet. Integral equation was solved numerically by the 
method based on the approximation of the integral using 
the quadrature formula with unequally spaced mesh. 

The comparison of the suggested model with our ear- 
lier model developed for gas absorption by large droplets 
(Kn <C 1) (see, Elperin eta al., 2008) showed that both 
models yield the same time of complete saturation of a 
large droplet by the soluble gas. 

Using the suggested model we studied absorption of sul- 
fur dioxide (SO2), dinitrogen trioxide (N2O3) and chlorine 
(CI2) by water nanoaerosol. It is showed that enhanced 
depiction of the dissolved N2O3 gas in a water droplet due 



to chemical reaction leads to the decrease of N2O3 concen- 
tration in the bulk of a water droplet and to the increase 
of the concentration gradient at the gas-liquid interface. 
Consequently, the flux of dinitrogen trioxide into a droplet 
is higher than the fluxes of sulfur dioxide and chlorine at 
later stages of gas absorption. 

It was demonstrated that neglecting kinetic effects leads 
to the significant overestimation of the soluble gas flux into 
a droplet during all the period of gas absorption. 

Appendix A. Derivation of the integral equation 
of Volterra type for the molecular 
flux density 

Here we give the details of derivation of Eq. (TT5j) . Let us 
first remove j from the boundary condition (|14[) . To this 
end let us introduce the new unknown function C(r,t): 



n L =C(r,t) 



Mr 

Dr ' 



(A.l) 



Substituting Eq. (|A.1[) into Eq. (ITT|) . and taking into 
account that in spherical coordinates the Laplacian Ar = 
2/r we obtain the following equation for C(r, t): 



dC_ 
~dt 



D L AC 



2j 



AC 



D L 



with the following boundary conditions 

dC 
dr 



and 



Substitution 



dC 
dr 



r=0 



C(r,t) 



J 

D L r=0 

X(r,t) 



reduces Eq. (|A.2|) to 



The boundary conditions to Eq. (IA.6|) read: 

= 



dx 
dr 



and 



Let us introduce the eigenfunctions 

d 2 (j) n 



dr 2 



= -K <p„ 



(A.2) 

(A.3) 

(A.4) 
(A.5) 

(A.6) 

(A.7) 
(A.8) 

(A.9) 



where the boundary conditions to the equation (|A.9|) are 
the same as for \, i-e., given by Eqs. (|A.7I) - (|A.8[) : 



d<j>n 

dr 



0, 



lr=0 



0. 



(A.10) 



Then the solution of the Eq. (|A.9|) reads: 

T 

4> n = u„sin(^„-) , 



(A.ll) 



where fi = na is the infinite set of the roots of the charac- 
teristic equation: 

/i = tan(/f) . (A. 12) 

The roots \i n of the Eq. (IA.12I) can be calculated numer- 
ically and are as follows /ii = 4.4934, = 7.7253, /13 — 
10.9041,^4 = 14.0662, = 17.2208 etc. The orthogonality 
condition for eigenfunctions reads: 



idr = 5 n 



(A.13) 



where S nm is the Kronecker delta. Equation (|A.13|) allows 
to determine the normalization constant u n : 

-1 



sin 2 (n n -)dr 



a sin n n 



(A.14) 



The eigenvalue fi n = (which is also the solution of Eq. 
(|A.12p ) and the respective eigenfunction require a special 
consideration. Solution of Eq. (|A.9p for //„ = reads: 



<p (r) = u Q r 



(A.15) 



where the normalization constant uq is determined from 

Eq. (TA~T3l) : 

3 



Uq 



(A.16) 



Let us now look for the solution to Eq. (IA.6I) in the fol- 
lowing form: 



X 



(r, t) = 2*o (t)4>o + ^n(t)Mr) , (A.17) 



n>0 



where the coefficient 2 appears due to double degeneration 
of the eigenvalue fj, n = 0. Equation for *„ reads: 

= -<r n * n + -^-(jt + Xj)b n - 2ja n , (A.18) 



where 



a n = X+D; 



a n = u n J sin 




in (/"n^J dr, 



sin 



dr 



(A.19) 



and 



(T = A, a = u I rdr = -V3a, 



bo = uq I r dr — -v3a 5 . 



(A.20) 



For n > we obtain 

Kj{t) ( /fl n ^ 2 



* n (i) 



and 



(^) &„+2 a „j y j(f) e -'-(*-*')df 

(A.21) 

2a / i(t')e~ A(t ~*'W (A.22) 



for n = 0. The number density of molecules A inside the 
particle near the surface reads: 



n_(a, i) 



x(M) j(*)o 



a D L 

Now let us determine \ from Eq. (|A.17|) . Noting that 



(A.23) 



b„<f> n (a) = 



n>0 



a 

T 



and 



f(^)\ + 2a„)<^(a) = 2 n>0, 
f(~) b n + 2a„ S j (f> n (a) = 3 n = . 



Eq. ((ATT7|> and Eq. (|A~23)) yield: 



(A.24) 

(A.25) 
(A.26) 

-^-^jit'jdt'. 



n>0 



(A.27) 

According to the Henry's law n + = Using the equa- 

tion for the flux density 

j(t) = -^a{a){ noo - n+) (A.28) 

we arrive at the following integral equation of Volterra type 
for j(t): 



m = 



a(a) 
Ana 2 



noo -HjS(t- t')j(t')dt' 




(A.29) 



where the kernel S(t — t') is given by the following formula: 
S(£) = 2 e-[^"/ Q ) 2+A l« + 3e- A « . (A.30) 



n> 
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